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Abstract. The accurate and efficient computation of the deformation of crystalline solids requires 
the coupling of atomistic models near lattice defects such as cracks and dislocations with coarse- 
grained models away from the defects. Quasicontinuum methods utilize a strain energy density 
derived from the Cauchy-Born rule for the coarse-grained model. 

Several quasicontinuum methods have been proposed to couple the atomistic model with the 
Cauchy-Born strain energy density. The quasi-nonlocal coupling method is easy to implement and 
achieves a reasonably accurate coupling for short range interactions. In this paper, we give a new for- 
mulation of the quasi- nonlocal method in one space dimension that allows its extension to arbitrary 
finite range interactions. We also give an analysis of the stability and accuracy of a linearization of 
our generalized quasi-nonlocal method that holds for strains up to lattice instabilities. 
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1. Introduction 

The motivation for multiscale coupling methods is that the accuracy of a fine scale model is often 
only needed in localized regions of the computational domain containing defects, but a coarse- 
grained continuum model is needed to make efficient the computation of large enough systems. 
The first quasicontinuum approximation was the energy-based quasicontinuum model (denoted 
QCE [15J). However, there are interfacial forces (called "ghost forces") when a uniform strain is 
modeled by the QCE energy [2JI18J . or equivalently, the uniform strain is not an equilibrium solution 
for the QCE energy (even though the uniform strain is an equilibrium for the purely atomistic and 
coarse-grained continuum models). The effect of the ghost force on the error in computing the 
deformation and the lattice stability by the QCE approximation has been analyzed in [2'U41I8],IT1]. 

Thus, there is a need for a more accurate atomistic-to-continuum coupling method. The ghost 
force correction method (GFC) achieves an increased accuracy by adding a correction to the ghost 
forces as a dead load during a quasistatic process [21 El [91 [13l [TH] - More accurate coupling can be 
achieved by a force-based approximation [ H I5 1 I7 1 [T2] , but the non-conservative equilibrium equations 
make the iterative solution and the determination of lattice stability more challenging |9j . 

An alternative approach is to develop a quasicontinuum energy that is more accurate than QCE. 
The quasi-nonlocal energy (QNL) was the first quasicontinuum energy without ghost forces in 
the atomistic-to-continuum interface for a uniformly strained lattice [19]. The developers of the 
QNL method introduced interfacial atoms, which interacted with the atomistic region using the 
atomistic model and interacted with the continuum region using the continuum model. For a one 
dimensional chain, the original QNL method is restricted to next-nearest-neighbor interactions. In 
this paper, we formulate the one-dimensional QNL energy in terms of interactions rather than the 
energy contributions of individual atoms [8] , which allows us to generalize the original QNL energy 
to arbitrary finite range interactions. A closely related method has been independently proposed 
and studied numerically for one and two dimensional problems by Shapeev in |17j . 
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We also give an analysis of the stability and accuracy of a linearization of the generalized quasi- 
nonlocal method in the one dimensional case. The stability and optimal order error analysis for a 
linearization of the original QNL model with second nearest-neighbor interaction was analyzed for 
the one dimensional case in [5j. A nonlinear a priori and a posteriori error analysis for the QNL 
model with second nearest-neighbor interaction in one dimension was given in |16j . 

In Section [21 we describe the coarse-grained local QC energy and the motivation for the devel- 
opment of more accurate QC energies. 

In Section O we present the notation used in this paper. We define the displacement space U and 
the deformation space 3^f- We then introduce the norms we will use to estimate the consistency 
error and the displacement gradient error. 

In Section HI we introduce the general QNL energy with finite range interaction and give an 
expression of the energy and equilibrium equations associated with the model. 

In Section [5j we give sharp stability estimates for both the fully atomistic model and the gen- 
eralized QNL model with respect to a uniform strain. Sharp stability estimates are necessary to 
determine whether quasicontinuum methods (or other coupling methods) are accurate near insta- 
bilities such as defect formation or crack propagation [6l[8] . Similar stability estimates of the fully 
atomistic and fully local quasi-continuum (QCL) models have been studied by discrete Fourier 
analysis in |llj . 

In section [6j we study the convergence rate of the generalized QNL model. We compare the 
equilibrium solution of the generalized QNL model with that of the fully atomistic model, and we 
use the negative norm estimation method [7j to obtain an optimal rate of convergence of the strain 
error. The error estimate depends only on the smoothness of the strain in the continuum region 
and holds near lattice instabilities, thus demonstrating that the generalized QNL method can give 
a small error if defects are captured in the atomistic region. 

2. The Local QC Approximation and its Error 

In this section, we first briefly describe the ID atomistic chain model. We give the fully atom- 
istic energy, £ a (Y), and then use linear interpolation to derive the coarse-grained local QC energy, 
£ qc {Y). We next compute and analyze the difference between these two energies. We then rescale 
the model and conclude that the coarse-grained local QC energy is formally a second order approx- 
imation (where the small parameter is lattice spacing scaled by the size of the macroscopic domain) 
to the fully atomistic energy when the strain gradient is small. 

For simplicity, we assume that the infinite atomistic chain with positions ye < ye+i has period 
2N, that is, 

ye+i+2N - ye+2N = ye+i -yi V£ e Z, 
where Z is the set of integers. We compute the energies of interaction using a potential <j) : (0, oo] — > 
R which we assume satisfies 

(1) <^GC7 4 ((0,oo],lR); 

(2) there exists r* so that <j)"(r) > for all r 6 (0,r*], and (f)"(r) < for all r £ (r*,oo); 

(3) (J) an d its derivatives decay at infinity. 

Examples of commonly used interaction potentials which satisfy the above conditions include 
the Lennard-Jones potential and the Morse potential. The total stored atomistic energy per period 
for the interaction potential, (j)(r), with up to sth nearest-neighbor interactions is given by 

s N 

£a (y)-=Yl Yl kvt+k-i-vi-i)- (2.i) 

k=ie=-N+i 
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We note that £ a (y) can be rewritten as 

s N 

£a (y) = Y Y HKve-ye-i)) 

k=l£=-N+l 

s N ( fc-1 ^| 

+ Y Y \ &(yt+k-i - vt-x) - r Y H k (y?+t - vt+t-i)) \ 

k=2l=-N+l { t=0 J 

N 

= Y & cb (ye ~ yz- 1 ) 

t=-N+l 

s N ( fe-1 "| 

+ Y Y \ ^(ye+k-i - vt-x) - t Y&i^yz+t ~ ye+t-i)) \ 



(2.2) 



where 



k 

k=2£=-N+l (. t=0 



W*0 := Y 

fe=l 



is the Cauchy-Born energy density [5|ll5j . 

Intuitively, we can reduce the total amount of computation by first choosing 2M, M « N, 
representative atoms in one period, linearly interpolating the positions of the remaining atoms, and 
then computing the total atomistic energy by using the interpolated positions. More precisely, we 
introduce the representative atoms with positions Yj such that 

Yj = y tj , j = -M,...,M, 
where the subindex £j for j = —M, . . . , M satisfies 

-N = 1-m < • • ■ < ij-i < £j < ■ ■ ■ < hi = N. 
The representative atoms thus satisfy 

Y-M <■■ < Yj-\ <Yj <■■■ < Y M . 

We denote the number of atoms between lj-i and Yj as Vj, and the distance separating the Vj 
equally spaced atoms between K'-i and Yj as fj, that is, 

Vj ;= £j - tj^X and f j = (Y j -Y j ^ 1 )/v j for j = -M + 1, . . . , M. 

Then the positions of the atoms between lj-i and Yj can be approximated by 

Vij-i+i = Vtj-to-t) = Y 3-i + irj = - 1 7 —Y j „ 1 + —Y j} 0<i<vj. (2.3) 

We thus have that 

Yj-i = ye J _ 1 < yij_ 1+ i <■■■ < y^-^-i) < Vtj = Y h < i < £/_,-. 

We further assume that each Vj is sufficiently large, i.e., Vj > s for j = —M + 1, . . . , M, so that 
the the span of a bond can cross at most one interface. We can then observe for the atomistic 
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deformations given by the interpolation (|2,3p that 
fe-1 

4>(ye+k-i - yi-i) - -r E HHvi+t - yi+t-i)) 

(o, £ j - 1 <£<£ j -(k-2), 

= \4>{(p + l)fj + (k-p- l)f i+ i ) - egfofj) - ^=^4>(kr j+1 ), i = £j-p,0<p<k-2. 

Rearranging the terms in (12, 2j) and applying the equalities above, we can rewrite the total atomistic 
energy £ a (Y) as the sum of a coarse-grained local QC energy, £ qc (Y), and an interfacial energy: 

s N 

£ a (Y) = J2 E HHyt-ve-i)) 

k=ie=-N+i 

s N ( fe-1 "J 

+ E E \ <K^+*-i - Ve-i) - r E 0( k (ye+t - Vl+t-i)) \ (2-4) 

fc=2^=-Af+l I i=0 J 

A/ 

= £ qc (Y) + £ ^> 

j=-M+l 

where the coarse-grained local QC model is 

M 
£ qc (Y):= Yl vjkbifj), (2-5) 

j=-M+l 

and where the interfacial energy is 

s fe— 1 , 

V i := E E {^ ^ f J + ( k " P)^'+ 1 ) " 1 4>( kf i) ~ -r^(*f>+i)}> J = -M + 1, . . . , M. 

fc=2p=l 

We note that we have used the periodicity for the representative atoms chosen in one period. 
Therefore, we have that fjvf+i = f_M+l and the definition of Vm makes sense. 

To understand the difference between the atomistic energy £ a (Y) and the coarse-grained local 
QC energy £ gc (Y), we first evaluate the interfacial energy terms, namely Vj, j = —M + 1, . . . , M. 
Using the Lagrange form of the Taylor expansion of degree 2 at |(fj + i + fj), we can obtain that 
Vj is given by 

V 3 = E\ 12 +fc ^( fc (Wj + 0- ~ Vjk)r j+ i)j(r j+ i - rjf\ , j = -M + 1, . . . , M, (2.6) 

and rjjk is a scalar satifying < rjjk < 1. 

2.1. Scaled Models. We next consider a scaled version of the atomistic and local QC energies. 
Thus, we define the two scaled interaction potentials 4>(r) and 4> c b{r) by 

4>(r) = -<j>(re) and 4> c b(r) = -4> cb (re) 
where e > scales the reference lattice. We can now convert (12.61) to the scaled form: 



V i = e El fc 12 + V (fc(Wj + (1 - Vjk)r j+ i)) (r 3 - r J+1 A (2.7) 
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where Tj := fj/e and rjjk is a scalar satifying < rjjk < 1- 

Therefore, the total atomistic energy (j2.4|) has the scaled form 

M M 

S a (Y)= j; *,■&*(*,) + £ ^ 

j=-A/+l j=-M+l 

= E v i M r d) + e £ ( r J - r i+ 1 ) 2 E \ — To — ^" ( fc fefe r J + ( X ~ Vjk)rj+i) J [ 

j=-M+l j=-M+l fc=2 ^ ^ ' ^ 

M 

= £^(Y)+ £ H j+1 Cj {eH j+1 (Y*' +1 f} 

j=-M+l 

(2.8) 



where 



and 






X,. : =^., Hj^iXj-X^fo 



Y : 



/._ y i"^-i 



v - v 

n 3 

Proposition 2.1. The coarse-grained local QC energy £ 9C (Y) is formally a second order approxi- 
mation to the fully atomistic energy £ a (Y), or more precisely, 

M 

f(Y) = £*(Y) + £ ^i+iC^ {cfli+i(l3+i) a } . 

i=-M+i 



where -ff = max./ flj . 



3. Notation 



Before further discussing the scaled models, we present the notation used in this paper. We 
define the reference lattice 

eZ := {e£ : £ £ Z}, 

where e > scales the reference atomic spacing and Z is the set of integers. We then deform the 
reference lattice eZ uniformly into the lattice 

FeZ := {Fe£ : £ G Z}, 

where -F > is the macroscopic deformation gradient, and we define the corresponding deformation 
Yf by 

(yF)e ■= Fe£ for - oo < £ < oo. 
For simplicity, we consider the space IA of 2iV-periodic zero mean displacements u = {ui)n^.% from 
y F given by 

N 
U := {u : Ui + 2N = U£ fori? G Z, and V, n ^ = 0}i 
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and we thus admit deformations y from the space 

yF '■= {y : y = Yf + u for some u £ U}. 

We set e = 1/N throughout so that the reference length of the periodic domain is fixed. 
We define the discrete differentiation operator, Du, on periodic displacements by 

{-Uu)£ : = , — oo < I < oo. 

We note that (-Du)* is also 2A r -periodic in £ and satisfies the mean zero condition. We will denote 
(Du)p by Du£. Then, we define 

L> w u := , — oo < £ < oo. 



e 

We can define (Z)( 3 'u)„ and (Z?^ 4 ^u)„ in a similar way. To make the formulas concise, we also 
denote Dug by u'p, D^ 2 'ui by u", etc., when there is no confusion in the expressions. 
For a displacement ueW, we define the discrete £ 2 norm by 



ll u lk := e Yl l n ^l 2 

\ e=-N+i J 
The inner product associated with the £ 2 norm is 

N 

(v,w) := e 2, v e w e f° r an v,w € W. 

4. The Generalized Quasi-Nonlocal Approximation. 

From (|2.8p , we find that the difference between the atomistic energy and the coarse-grained local 
QC energy is formally of second order, i.e., maxj O(eHj). However, when the strain gradient is 
large in some regions, namely where eHj(Y!') 2 is large, this error can be unacceptable. Hence, to 
maintain both accuracy and efficiency, we should use the atomistic model where the strain gradient 
is large and the local QC model where the strain gradient is moderate. In this section, we propose 
a hybrid atomistic-continuum coupling model that extends the quasi-nonlocal model [19] to include 
finite-range interactions. To simplify our analysis, we will propose and study our quasi-nonlocal 
model without coarsening, although the local QC region can be coarse-grained. 

According to our assumption in Section [21 the total stored energy of the atomistic model per 
period (that includes up to the sth nearest neighbor pair interactions) in dimensionless form is 

N s /fc-1 \ 

£=-N+l k=l \j=0 J 

Here, the atomistic energy is a sum over the contributions from each bond and we can also rewrite 
it in terms of energy contributions of each atom, 

N 

£ a (y)=e Y, E e(y), where 
e=-N+i 

eky) := ^E^E^) + ^E^E^)- 

fe=l j=l fc=l j=0 
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If y is "smooth" near yg, i.e., y' i+ - and y' e _- vary slowly near y^ then we can accurately ap- 
proximate the distance between A;th nearest-neighbors of yi by that of first nearest-neighbors to 
approximate E^(y) by E%(y), where 

E e(y) ■= \ E Wti + </>M+i)] = \ [<p cb {y'i) + <t>Mi + i)] ■ 
fe=i 

If y'^ is "smooth" outside of a region {— K, . . . , K }, where 1 < K < iV, then the original 
quasicontinuum energy (denoted QCE) [15] uses the atomistic energy Ef in the atomistic region 
A := {— K, . . . , K} and the local QC energy E9 in the continuum region C := {— N + 1, . . . , N} \ A 
to obtain 

-A'+l K N 

£ qce (y)--=e E E ?(y) + e E ^ a (y) + e E ^ c (y)- 

£=-iV+l ^=-iC ^=iC+l 

Although the idea of the QCE method is simple and appealing, there are interfacial forces (called 
ghost forces) even for uniform strain 0H2]- The subsequent low order of consistency for the QCE 
method has been analyzed in jUEKll]. 

The first quasicontinuum energy without a ghost force, the quasi-nonlocal energy (QNL), was 
proposed in [19] . The QNL energy introduced in [19] was restricted to next-nearest neighbor 
interactions which in ID are next-nearest interactions. By understanding this method in terms of 
interactions rather than the energy contributions of "quasi-nonlocal" atoms, we have been able to 
extend the QNL energy (for pair potentials) to arbitrarily finite-range interactions while maintaining 
the uniform strain as an equilibrium (that is, there are no ghost forces). A generalization of the 
QNL energy to finite-range interactions from the point view of atoms was proposed in [10], but 
the construction requires the solution of large systems of linear equations that have so far not 
permitted a feasible general implementation. The interaction-based approach that we give here has 
been generalized to two space dimensions in |17j . 

In our QNL energy, the nearest-neighbor interactions are left unchanged. A A;th-nearest neighbor 

interaction (f> ( J2j=o Ve+j ) where k > 2 is left unchanged if at least one of the atoms £ + j : j = 
— 1, . . . , k — 1 belong to the atomistic region and is replaced by a Cauchy-Born approximation 

fe— l \ fe— l 

K j=0 J 3=0 

if all atoms £ + j : j = — 1, . . . , k — 1 belong to the continuum region. We define 
A qn i(k) := {-K -k + l,-K-k + 2,...,K,K + l}, 

C qn i(k) :={-N + l,...,-K-k}\J{K + 2,...,N} 

for k = 2, . . . , s. Then the our generalized QNL energy is given by 

TV s / k—1 \ s fe— 1 

£^(y):=e £ *(«$) + e£ E * £«ki + c £ E lE^W" ^ 
e=-N+i k=2£ £ A qnl (k) \j=o J k=2iec qnl (k) j=o 



(4.3) 



Since the forces at the atoms for the Cauchy-Born approximation (|4.2p are unchanged for uniform 
strain, the QNL energy does not have ghost forces. 

Proposition 4.1. The QNL energy defined in (|4.4p is consistent under a uniform deformation, 
i.e., it does not have a ghost force. 
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is consistent under a uniform deformation. We note that 

fc-l 



Proof. The force at atom £ is defined to be — ac \ (y> , so from (|4.4p . we only need to verify that the 
approximation 

8 



dyi+r 



'k-i 

J=0 







y=yF 



dye+m 



lEHkyl 



3=0 



Vm 



■l,...,k. 



y=yF 



Therefore, the QNL energy is consistent. □ 

We note that the negative of the QNL forces for a general deformation y qnl £ y F is given by 

N s /fc-l \ /fc-l 



{6SV*(y**),w) = e £ </>>(y'M + e £ £ 0' £ ^ £ «£tf 



€=— JV+l 



fc=2 £eA„ ni (fc) \t=0 



vi=0 



fc-l 



(4.5) 



+ e £ £ E^'M+jK+i VwG ^- 

fc=2 teC gn! (fc)j=0 

5. Sharp Stability Analyses of Atomistic and QNL model. 

A sharp stability analysis for both the atomistic and QNL models are needed to determine 
whether the QNL approximation is accurate for strains near the limits of lattice stability. In this 
section, we will thus give stability analyses for both models. In order to provide clear statements 
and proofs, we first apply a similar method given in [8] to a third nearest-neighbor interaction range 
problem and then generalize the results to the finite range case. 

5.1. Atomistic Model with Third Nearest-Neighbor Interaction Range. The total energy 
for the atomistic model given by a third nearest-neighbor interaction without external forces is 

N 

£ a (y) = £ £ M) + M + y , e+i) + M + yi+i + y'e + 2)- (5-i) 

e=-N+i 

It is easy to see that the uniform deformation y^v is an equilibrium of the atomistic model, that is, 

(5£ a (y F ),w) =0 VweW. 

We will say that the equilibrium y F is stable for the atomistic model if 5 2 £ a (y F ) is positive definite, 
that is, if 



N 



(5 2 £ a (y F )u,u)=e £ {4>' F \u' t \ 2 + <& F \u' t + u' e+1 \ 2 

e=-N+i 

+ <l>3 F \u' e + u' e+1 + u' e+2 \ 2 ) >0 VueW\{0} 



(5.2) 



where 

We observe that 



<t>" F : = 0"(n <t>2F ■= <t>"vn ■ ■ -a" f ■= <p'\sf). 



i,' -7/ 1 2 

u e+i u i\ 1 



\u' t + u' t+l \ 2 = 2\u',\ 2 + 2\u', +l \ 2 
K + u' e+1 + u' e+2 \ 2 = 3\u' e \ 2 + 3|*4 +1 | 2 + 3\u' e+2 \ 2 



(5.3) 



3 e 2 |< +1 | 2 - 36 2 |< +2 | 2 + e 2 K +2 - < +1 | 2 . 
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Because of the periodicity of u'g, u'g in £, we note that the sum of u' i+l from — N + 1 to N is equal 
to the sum of u' e and the sum of u" +1 is equal to the sum of u", etc. So we can rewrite (|5.2p and 
get the lower bound 

N N 

(5 2 £ a (y F )u,u)=e £ (^ + A^ F + 9^ F )\u' e \ 2 - e E (e 2 $ F + 6e 2 <^)K| 2 

e=-N+i e=-N+i 

N 

+ e E ( e 2 ^V)K + i-<| 2 

^=-Af+l 

(jjl \ ajJI i ni» \ II 'I|2 / 1 ±11 , o 2 j// \ ||„// 1| 2 

= (<Pf + 4 <t>2F + WsF) ll u 11*2 - (e <f>2F + 2e <P3f) \\ u \\% 

N 

-<^ F ) E [4KI 2 - K+i - <l 2 ] 

v. /j.// i ji» i ni" \ 11../ II 2 2/i// i oj." \ll "II 2 

> (<Pf + 4( P2F + y<P3F) ll u 11*2 - e (</> 2F + 20 3F J ||U || |2 



(5.4) 



AT AT AT 



and since 

i 
E 4|<| 2 > E [4|<| 2 -|< +1 -<|V E [2|<| 2 + 2< +1 <]>0, (5.5) 

e=-N+l £=-N+l l=-N+l 

the upper bound of (5 2 £ a (y f)u, u) is 

(5 2 £ a (y F )u, u) < ($, + 44% F + 9<% F ) ||u'||| - e 2 (^ F + 6^ F ) ||u"|||. (5.6) 

We thus define 

4 3> _ tf + 4& + 9& , „ := ^ »£k _ « , ( , 7) 

where the equality /i e = sin l 7re / J j s given in [8]. We then obtain the following stability result for 
the atomistic model. 

Theorem 5.1. Suppose 2 ' F < 0, <f>'^ F < and (f)'£ F = /or fc > 4. Then yp is a stable equilibrium 
of the atomistic model if and only if A F —e 2 fi 2 [<f>2 F +riF(t>^ F ) > 0, where T)f is a constant 2 < tjf < 6 
i/toi is uniquely determined by the deformation gradient F, and A F and fi e are defined in (|5.7p . 

Remark 5.1. The hypotheses (ft^p < 0, 3 '^ < in Theorem 15.11 are reasonable according to the 
assumptions of the interaction potential 4> m the Section [2J We note that y' e < 4- only occurs under 
extreme compression, and in that case the finite range pair interaction model (|5.ip itself can be 
expected to be invalid. We refer to Section 2.3 of [16J for further discussion of this point. 

5.2. The QNL Model with Third Nearest-Neighbor Interaction Range. The total energy 
for the QNL model given by a third nearest-neighbor interaction model without external forces is 



s qnl (y) = e e <^) + e E +M + M+* E liHiye) + HH + i)} 

£=-N+l £&A qnl (2) teC, ni (2) 



E <l>(vi + i/t+i + i/t+2)+e E oi^y^ + HSy'i+J + HSy'e^)}- 



(5.8) 



teA qnl (3) £&C qnl (S) 
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Since the QNL model does not have ghost forces, the uniform deformation y^ is still an equilibrium 
for (|5.8p . We thus focus on £™;(y.F)[ u > u ] to analyze the stability of the QNL model. We note that 
£'' nl (y F )[u, u] can be written as 

JV 
(5 2 £^ l (y F )u,u)=e £ <j>' F \u' e \ 2 

e=-N+i 

Or- A in\ or-r* ln\ v - * 



eeA qnl (2) tec qnl (2) 

+ e X] $>f K + <4+i + «/+2 1 2 

^€A,„,(3) 

Applying (|5.3p to f)5.9|) . we obtain 

(6 2 £« nl (y F )u,u) = 4 3) ||u'||| - e 3 £ ^VK+il 2 



(5.9) 



teA,j ni (2) 

e 3 £ ^{K + il 2 + K +2 l 2 + K +2 + ^ + il 2 } (5-10) 

(3) 1. /,,2 



>A^||u'| U? . 

Since the lower bound in (|5.10p is achieved by any displacement supported in the local region (which 
exists unless K €. {N — 2, N — 1, N}), it follows that y F is stable in the generalized QNL model if 

and only if A$ > 0. 

Theorem 5.2. Suppose 2 V — 0> ^V < and ^^ = for k > 4. XTten y_p is a stable equilibrium 

(3) 

0/ i/ie QNL model if and only if A F > 0. 

We have given above a stability analyses of the atomistic model and the QNL model with a 
third-nearest neighbor interaction range. We now study the general case. 

5.3. The Atomistic Model with sth Nearest-Neighbor Interaction Range. In the case of 
sth nearest neighbor interactions, the total energy for the atomistic model without external forces 

is 

N s /fe-1 \ 

£a (y) = e £ X> £«ki ■ (5-n) 

e=-N+i fc=i \j=o J 
The uniform deformation y^ is still its equilibrium, so the stability condition is 

N s /fe-1 \ 2 

(S 2 £ a (y F )u,u) = e £ E^ E%i >0 Vu € W\ {0}. (5.12) 

£=-N+lk=l \j=o I 
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We assume that (f>" F > and (j)'£ F < for k = 2, . . . , s. We first consider (f)^ F [Ylj=o u i+i ) ^ or 



k> 2: 



fc-1 \ fc-1 fc-2 fe-1 

b kF ( E U 't+J = ^kF Y K+i) 2 + 0'kF Y Y 2u 't+3 U 't+i 
,j=0 / j=0 j=0i=j+l 



fe-1 



fc-2 fc-1 



b 'kF Y ( u 't+j) 2 + ^feF Y Y { ( u 'e+j) 2 + K+if - K-h - u i+i) 2 } 



i=o 



j=0 j=j+l 



(5.13) 



Recalling that u'l :- 



'fc-i 
^fcF ( 2^/ u t+j 

J=0 



— , we can further simplify (|5.13p and get 

2 



fc-1 



fc-2 fc-1 



^Efef+^E E {K+/ + Kh) 2 } 



i=o 



i=0 i=j+l 



fc-2 fc-1 / i 

-^ 2 EE E< + * 

i=o<=i+i V=j+i / ( 514 ) 

>^E(« 8 +^E E {(^) 2 + K +i ) 2 } 

i=0 j'=0 i=j+l 

fc-2 

-^F^^K+j+i) 2 - 

The last inequality comes from the fact that — (ft'lpt 2 > for /c = 2, . . . , s, so we determine that 
the terms i = j + 2, . . . , k — 1 in the third expression above are all nonnegative. Therefore, (j5.12[) 
becomes 

N 

(5 2 £ a (y F )u,u)>e Y <t> F { u ef 

fc-AT+l 



s AT 



fc-1 



fc-2 fc-1 



+ E e E ^|E(i/+EEN 2 +(4) : 

Y k2 ^F\\A%-Y^F(k-l)\\u" 



fc=l 



j=0 j=0 t=j+] 

^fcM*-w"" 2 

fc=2 



fc-2 

^EK +J - 

j=0 



(5.15) 



On the other hand, we can further study the third term in the fist line of (I5.14p and rewrite it as: 

2 r 

fc-2 fc-1 

E u 'i+t 2 + Y Y ( u i+t 2 + u e+s 2 ) 

t=j+l t=j+l s=t+l 

+^ 2 EE E E W +t -uW 

j=0 i=j+l t=j+l s=t+l 



fc-2 fc-1 

^Y Y | E <* 

j=0i=j+l \t=j+l 



fc-2 fc-1 

^Y Y 

3=0 t=i+l 
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Thus, we can use the first line of (|5,14p and (f>'L F < 0, k = 2, . . . , s to obtain an upper bound of 
0fcF (E*=o *4+j) : 

( k— 1 \ k— 1 A;— 2 k— 1 

^ E^+i <&EK + /+^EE {K+/ + K+*) 2 } 



J=o 



j=0 j=0 i=j+l 

k-2 k-1 

^E E 



j=0 i=j+l 



i i— 1 i 

E n "+* 2 + E E ( u "+t 2 + u 'U 

t=j+i t=i+i s=t+i 



(5.16) 



Observing the last term of (|5.16p and because of the periodicity condition of u", for each 



^fc-2 v^fc-l 



fixed k, we have that £ i=0 E»=j+i Et=j+i n "+* + Et=,'+i El=t+i K+ 1 + n : 



l ^+t ^ "£fs 



is equiva- 



^fc-2 v^fc-i 



lent to £ i=0 Ei=j+i(* - i) l u 

C(yF)[u,u]: 



2|„,//|2 _ fc 4 -fc 2 i ,//|2 



12 



u»| . Thus, we can obtain the following upper bound of 



(5 2 £«(y F )u,u) <E fc V^iKii| -E e V,vE E (*-^) 2 E e i<i 

j=0 i=i+l £=-7V+l 



fc=i 



fc=2 



E/2j// II /||2 V^ 2 ill 



k=l 



Note that when k = 2, we have — ^ — 



fc=2 



fe 4 - k 2 
12 



(5.17) 



,"l|2 



1 = A; — 1, which means that the upper bound equals the 



lower bound if s = 2. 
We next define 



AJ. :=£*%., 



(5.18) 



fc=i 



and use the same notation for /i e defined in (|5.7p . We then obtain the following sharp stability 
estimate: 

Theorem 5.3. Suppose (f>'l F < for k = 2, . . . , s. There exists a constant B = Bp satisfying 

£> - m F >b f > <&, + J2 ^r^ 

k=2 fc=3 

such that yp is a stable equilibrium of the atomistic model if and only if A F — e 2 fi 2 Bp > 0, where 
A S F is defined in (J5.18P and fj, e is defined in ([57 



Remark 5.2. The assumptions (j>'£ F < for k = 2,...,s in Theorem 15.31 are reasonable by the 
considerations given in Remark 15.11 

5.4. The QNL Model with sth Nearest-Neighbor Interaction Range. We now consider 
the stability of the QNL model in the general case. Again y^ is an equilibrium of the QNL model 
(|4.5p when there is no external force, so we need to check whether 



(8 2 £i nl (y F )u,u)>0, VuGW\{0}, 
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where {5 2 6 qnl (y f)u-, u) is 

(5^\y F )u,u)=e £ <j>' F {u' E f + e^ E ^ ( X>*fi 

«=-iV+l fc=2£ S A 9ni (fc) \j=0 

s k—1 

+ e E E ^^K +i ) 2 

k=2e & C qnl (k) j=0 

N 



e E <^04) 5 



(5.19) 



-JV+l 



fc-1 



fc-2 fc-1 



+ e E E tf* E K+/ + E E [K+i) 2 + «h) 

fc=2te^ gni (fc) [j=° j=0i=j+l 

(\ 2\ 
i \ s fc-1 

E < + « 
t=j+i 



> + e E E fc E^K+i) ' 

k=2ieC qnl (k) j=o 



and for k > 2 



A qn i(k) = {-K -k + l,-K-k + 2,...,K,K + l}, 

C qn i(k) = {-N+l,...,-K-k}\J{K + 2,...,N}. 

Since (f/£ F < when k > 2, ()5.19j) becomes 

AT s r A:— 1 

(j 2 ^(y F ) U ,u) = e e ^K) 2 + ^E E ^EW J 

«=-iV+l fc=2£ Gj 4 gni (fc) li=0 



fc-2 fc-1 



fc-2 fc-1 



+££ [(«l t ,f+W']-EE<' E 

j=0i=j+l j=0i=j+l \t=j+l 

+ ^E E ^E^K +J ) 2 

k=2£eC qnl (k) j=0 

= AU\u'\\%-e±e^ F E EE^E^ 

fc=2 £ Gj4gn! (fc)j=0i=i+l \t=j+l 

> 4 s lli/ll 2 

— -^-fII" Hi 12 ' 



u i+t 



(5.20) 



where A^ is defined in (|5.18p . Since the lower bound in (|5.2U|) is achieved by any displacement 
supported in the local region (which exists unless K G {N — s + 1, . . . , N}), it follows that y F is 
stable in the QNL model if and only if A S F > 0. 

Theorem 5.4. Suppose that K < N — s + 1 and that (f)'l F < for k = 2, . . . , s. Then y F is a stable 
equilibrium of the QNL model if and only if A F > 0. 
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Remark 5.3. From Theorem 15.31 and Theorem 15.41 we conclude that the difference between the 
sharp stability conditions of the fully atomistic and the QNL models is of order 0(e 2 ). This result 
is the same as for the pair potential case in [6j. 

6. Convergence of the QNL model. 

So far, we have investigated the stability of the fully atomistic model and the QNL model. In 
this section, we will give an optimal order error analysis for the QNL model. We compare the QNL 
solution with the atomistic solution and give an error estimate in terms of the deformation in the 
continuum region. 

6.1. The Atomistic Model with External Dead Load. The total atomistic energy with an 
external dead load f is 

£? ot {y)-.= £ a {y)+Hy) Vye^, 

where 

N 

F(y) : = - E e ^ ye -- 

e=-N+i 
To guarantee the existence of energy-minimizing deformations, we assume that the external loading 
force f is in ti. The equilibrium solution y a G yp of the atomistic model with external force f then 
satisfies 

-(5£ a (y a ) lW ) = (5T(y a ),w) VweW, (6.1) 

where 

TV s k—1 /k—1 \ 

£=-N+lk=lj=0 \t=0 J 

and the external force is given by 

N dT N 

(&F(y),w) := ^J —(y) We = - ^ efgWf 

e=-N+i yi e=-N+i 

6.2. The General QNL Model with External Dead Load. The total energy of the QNL 
model corresponding to a deformation y 6 yp is 

£?:l(y)--=£ qnl {y)+Hy) 

So, the equilibrium solution y qnl S 3^7 satisfies 

- {5£ qnl (y qnl ),w) = {5T(y qnl ),w) Vw e U, (6.3) 

where 

N s /fc-1 \ /fc-1 

(6e**(y**), w) = 6 y, MM + * E E <M E y't+t E «&« 

e=-N+i k=2 ieA qnl (k) \t=o J \j=o , , 

s k—1 

+ e E E E^M+^K+j VwGW - 

fc=2 i£C qnl (k)j=0 

Setting y qnl = yp + u qnl and y a = y^ + u a , where both u qnl and u a belong to U, we define the 
quasicontinuum error to be 

e qnl ._ a _ qnl _ ^a _ ^qnl ^ 
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To simplify the error analysis, we consider the linearization of the atomistic equilibrium equations 
(|6.ip and the associated QNL equilibrium equations (J6.3P about the uniform deformation y F . The 
linearized atomistic equation is 

- (5 2 £ a (y F ) u a , w) = (ST(y F ),w) for all weW, (6.5) 

and the linearized QNL equation is 

- (5 2 £ qnl (y F ) u qnl ,w) = {5F(y F ),w) for all w£«. (6.6) 

We thus analyze the linearized error equation 

(5 2 £ qnl (y F ) e qnl ,w) = (T^, w) for all weU, (6.7) 

where the linearized consistency error is given by 

(T^, w) := (5 2 £ qnl (y F ) u a , w) - (5 2 £ a (y F ) u a , w). (6.8) 

To obtain an optimal order consistency error estimate, we extend the negative norm method 
given in 0[16] to the sth nearest-neighbor linearized consistency error functional 

s fc-1 / fc-1 \ 

<T^,w> = e]T E <t>'l F Y.[ kDu ^-Y. Du z+A Dw z+v VweW. (6.9) 

k=2iec qnl (k) j=o V t=o / 

For each fixed k, we define T q k n := e^' F Efec nl (k) Sj=o [kDu i+ j - ^1$ Dui +t ). Then we 
have 

-K-k fc-1 / fc-1 \ 

(TrU) = C E E (kDu e+J -Y,Du e+t \Dw t+3 

N fc-1 / fc-1 \ 

+ e ^V E E kDu e+j - E Du ^ Dwi +i ( 6 - 10 ) 

£=K+2j=0 V 4=0 / 

2N~K-kk-l I fc-1 \ 

= ^If E E kDu e+j ~ E Dui + t Dw t+r 
£=K+2 j=0 V t=o / 

We use the 2N periodicity of the model in the last equality to simplify the expression of (|6.10p . 

Note that we can change the indices, rearrange the order of summation of the sums in the last 
equality of (|6.10j) . and rewrite it as the following expression 

2N-K-k k 

(Tf, w) = e<j% F E E( fc - ft i- Du e-i + 2Du t - Du i+i) Dw e 

i=K+k+l j=l 

K+k £-{K+2) / fc-l \ 

+ e( t > kF E E [kDu e -Y,Du e - j+t \Dw £ (6.11) 

e=K+2 i=o V t=o / 

2N-K-k fc-1 / fc-1 \ 

+ e( PkF E E kDut+j - E Du £+t Dwi+j- 

t=2N-K-2k+2j=2N-K-k+l~£ V t=0 / 

The first term of (|6.1ip corresponds to the inner continuum region, which is of second-order because 
of the symmetries of the interaction. The second and the third terms are the interfacial terms. They 
are only of first-order since they lose the symmetries of the interaction. 
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Now we will give an estimate of the consistency error (T 9 , w) in the following theorem. We 
first define the following semi-norms: 



v||;,„ ,.,:=€ Yl v l and ll v ll|(V/- ; - ( 



|2 



,l(k)) - 



teC 9n ((fc) 



E 

e&Z q „l{k) 



V$, 



where I qn i(k) := {—K — k + 1, . . . , — K — 1} |J{-ftT + 2, . . . , K + k} is the interface between the 
continuum and atomistic regions, C qn i(k) := C qn i(k) [jl qn i(k). 

Theorem 6.1. The consistency error (T qnl ,w), given in (|6.9p . satisfies the following negative 
norm, estimate 



(T^,w) < £e 2 Ci(AMiHP (3) u| 



^(C,nl(*)) 



,fc=2 



+ J2eC 2 (kM F \ (V2ii) ||D (3) u||^ (Vi(fc)) 1 ||Dw 



(6.12) 



1^; 



fc=2 



where Ci(k), C2(k) are positive constants independent of e. 
Proof. From (|6.1ip . we have 

(Tf ,w)| < e 2 C 1 (^)|^Vlll^ ,(3) u|U ?( c g „ i (A ; ))ll^ w ll^ +eC 2 (^)|^V|||,C»^) u | U?(x ^ (fc)) || jDw ||^. (6 . 13) 
Therefore, we obtain an optimal order estimate for (|6.9p 



Irr^qnl 



w) 



< 



^6 2 Cl (fc)|^ F IP (3) u| 



tl(C qn i(k)) 



k=2 



k=2 



We note that we have 



ll D,2, "ll|(V« (t » = « E ( D(2, «') s £ ||o (2 »u| 

te/(fc) 

<e2 S || J D( 2 )u||2 ?0(V;(A;)) . 



</4Vl ll^ (2) u||^(x,„Kfc)) 


ll-Dw 11/2. 

ii nt e 




(6.14) 


«S°(2^(ft)) 2^ X 





«,„,(/:) 



Thus, we can obtain from (|6.14p the more concise (but weaker) estimate 
(T^,w)| < \jZ^mK F \\\D^Ap A e qnim 

+ ^eC 2 (kM F \ (V2sl) ||Z^u|| Wgn;W) 1 \\Dw\\ 



.k=2 



k=2 



(6.15) 



□ 



We can finally derive the convergence rate of the generalized QNL model from the above consis- 
tency error estimate and the stability estimate (|5.20p . 



Theorem 6.2. Suppose that A F > 0, where A S F is defined in (|5.18p . Then the linearized atomistic 
problem (|6.5p as well as the linearized QNL approximation (|6.6p have unique solutions, and they 
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satisfy the error estimate 

\\Du* Dn«*\\ < g^!!^ fc )l ^ ll|g(3)u| ^(^W) 



4 s 



Zi=2e 3/2 (V2^C 2 (kM F \\\D®u\\e° 



(Ignl(k)) 



A F 

Proof. This error estimate for the generalized QNL model follows from the error equation (|6.7p . 
the stability result (|5.20[) . and the estimate of the consistency error (|6.15|> . D 

Remark 6.1. The above 0(e 3 ' 2 ) error estimate for the strain in £ 2 is the optimal order of conver- 
gence. For the next nearest neighbor interaction, an explicit expression of the linearized error is 
given in [5] that justifies this claim. 

Remark 6.2. We note that the error estimate in Theorem 16.21 depends only on the smoothness 
of the strain in the continuum and interfacial regions and it holds for linearizations near lattice 
instabilities according to the stability estimate in Theorem 15.41 thus the generalized QNL method 
can give a small error if defects are captured in the atomistic region. 

7. Conclusion. 

We propose a generalization of the one-dimensional QNL method to allow for arbitrary finite 
range interactions. We study the stability and convergence of a linearization of the generalized 
QNL energy with arbitrary sth nearest-neighbor interaction range. We extend the methods given 
in [6j[8] to give sharp conditions under which the atomistic model and the generalized QNL model 
are stable. The difference of the stability conditions between the QNL and atomistic model is 
shown to be of order 0(e 2 ). 

We then give a negative norm estimate for the consistency error and generalize the conclusions 
in [5] to the finite-range interaction case. We compare the equilibria of the generalized QNL model 
and the atomistic model and give an optimal order 0(e 3 ' 2 ) error estimate in £ 2 for the strain that 
depends only on the deformation in the continuum region. 
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